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Recently, a relativistic gravitation theory has been proposed [J. D. Bekenstein, Phys. Rev. D 70, 
083509 (2004)] that gives the Modified Newtonian Dynamics (or MOND) in the weak acceleration 
regime. The theory is based on three dynamic gravitational fields and succeeds in explaining a large 
part of extragalactic and gravitational lensing phenomenology without invoking dark matter. In this 
work we consider the strong gravity regime of TeVeS. We study spherically symmetric, static and 
vacuum spacetimes relevant for a non-rotating black hole or the exterior of a star. Two branches of 
solutions are identified: in the first the vector field is aligned with the time direction while in the 
second the vector field has a non- vanishing radial component. We show that in the first branch of 
solutions the /3 and 7 PPN coefficients in TeVeS are identical to these of general relativity (GR) while 
1^-^ , in the second the f3 PPN coefficient differs from unity violating observational determinations of it 

f***^ ■ (for the choice of the free function F of the theory made in Bekenstein's paper). For the first branch 

of solutions, we derive analytic expressions for the physical metric and discuss their implications. 
£\J , Applying these solutions to the case of black holes, it is shown that they violate causality (since 

, | they allow for superluminal propagation of metric, vector and scalar waves) in the vicinity of the 

event horizon and/or that they are characterized by negative energy density carried by the fields. 

^ . PACS numbers: 95.35.+d, 04.80.Cc, 04.70.-s 

(N 

I. INTRODUCTION 

(N 

^vq ■ On cosmological scales, Newtonian gravitational theory under-predicts the acceleration of stars and gas. Further- 
more, galaxies and clusters of galaxies show anomalously large gravitational lensing when only baryonic matter is 
taken into account. A natural "cure" for these discrepancies is to assume the existence of dark matter, which domi- 
nates over visible matter [H H> El S. § - Such dark matter might solve the "missing mass" problem within the standard 
theory of gravity (i.e., general relativity or GR). For this picture to be complete, however, the origin of dark matter 
needs to be identified. 

A second approach to the acceleration discrepancy and lensing anomaly is to look for alternative theories ofgravity, 
that modify GR on large scales. Among many other attempts, the MOND paradigm has been proposed [g,Ll- It is 
I characterized by an acceleration scale do so that 
5h ' 

jx(\a\/a )a = -V$n, (1) 



in 
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where $n is to be understood as the Newtonian potential, jl(x) — x for i<1, while fj,(x) — 1 for x ^> 1. 

This empirical law has been very successful in explaining the rotation curves in a large number of (spiral, low 
surface brightness and elliptical) galaxies using the observed distributions of gas and stars as input 0, 0, Inj. MOND 
can also explain the observed correlation between the infrared luminosity of a disk galaxy Lk and the asymptotic 
rotational velocity v a (i.e. the the Tully-Fisher law: Lk oc ), on the (suggested by population synthesis models) 

assumption that the M/Lk ratio is constant. 

However, MOND is merely a prescription for gravity and not a self consistent theory. It violates, for example, 
conservation of momentum and angular momentum and does not provide the formulation to describe light deflection 
or to build a cosmological model. A theory of gravity is needed that has the MOND characteristics in the weak 
acceleration limit but also has full predictive power. 

In this context, a new relativistic theory of gravity has been introduced by Bekenstein [j^. It consists of three 
dynamical gravitation fields: a tensor field {gap), a vector field (U a ) and a scalar field (<^>) leading to the acronym 
TeVeS. The theory involves a free function F, a length scale i (that can be related to ao) and two positive dimensionless 
constants k, K. 

TeVeS gives MOND in the weak acceleration limit (and therefore inherits the successes of MOND on large scale), 
makes similar prediction on gravitational lensing as GR (with dark matter) and provides a formulation for constructing 
cosmological models. As a drawback however, one can mention that Te VeS still appears to need dark matter to address 
the cosmological matter problem, i.e., the fact that observations require that the source term of Friedmann's equations 
is a factor of ~ 6 the baryonic matter density. 

A self consistent theory must be causal (i.e. not to allow for superluminal propagation of any measurable field or 
energy) and in TeVeS this is the case provided that <fi > ( Section VIII). It can be shown that, for a range of 
initial conditions, FRW cosmological models with flat spaces in TeVeS expand for ever with < cj) -C I throughout. 
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Moreover, in the vicinity of a star embedded in this cosmological background <f> is still positive. So, in a wide range 
of environments, TeVeS has been shown to be causal. 

The predictions of the theory have thus been explored to some extent for a range of strengths of the gravitational 
field: from the MOND limit to the Post-Newtonian corrections in the inner solar system. The strong gravity regime 
(e.g. in the vicinity of a black hole or a neutron star) of TeVeS has not been studied. This regime is a topic of this 
work. The motivation for this work is two-fold. First, it is interesting to see how the physics of compact objects (i.e. 
black holes, neutron stars) differ in TeVeS with respect to GR and what constraints (if any) observations can put 
on its free parameters. Second, one can check the consistency of the theory (e.g., its causality, positivity of energy 
carried by the fields) in these extreme conditions. 

In Section II, we summarize the fundamentals of TeVeS and in Section III wc consider its strong gravity limit. 
We limit ourselves to static, spherically symmetric and vacuum spacetimes relevant for a non-rotating black hole or 
the exterior of a star. Two branches of solutions are identified: in the first the vector field is aligned with the time 
direction while in the (not previously explored) second branch the vector field has a non- vanishing radial component. 
We show that the (3 and 7 PPN coefficients in TeVeS are identical to these of GR in the first branch of solutions 
while the (3 PPN coefficient differs in the second. For the choice of the free function F made in Ref. 0], we find 
that TeVeS predicts a value for (3 that is in conflict with recent observational determinations of it. In Section IV, we 
consider the first branch of solutions and derive exact solutions for the metric for arbitrary values of the parameters 
of the theory. The observational properties of the black holes in TeVeS are discussed in Sec. V along with the issue 
of superluminal propagation of waves in the black hole vicinity. Conclusions are given in Sec. VI. 

II. THE BASIC EQUATIONS OF TEVES 

TeVeS is based on three dynamical gravitational fields: a tensor field (the Einstein metric g a p), a 4-vector field U a 
and a scalar field <f> with an additional nondynamical scalar field a. The physical metric g a/ 3 in TeVeS is connected 
to these fields through the expression 

g aP = e- 2 ^g a0 - 2U a U sinh(20). (2) 

The total action in TeVeS is the sum of 4 terms S g , S s , Su and S m (see Ref. H3)) where S g is identical to the 
Hilbert-Einstein action and is the part that corresponds to the tensor field, while S s , Sjj, S m are the actions of the 
two scalar fields, the vector field and the matter respectively. The basic equations of TeVeS are derived by varying 
the total action S with respect to g a ^ , <j>, a, U a . 

Doing so for g a @ , one arrives to the metric equations 

G af3 = 87rG[f al3 + (1 - e^)U% {a U p) + r a0 ] + G a0 , (3) 

where a pair of indices surrounded by parenthesis stands for symmetrization, i.e. Ai a Bm = A a Bp + ApB ai the G a p 
denotes the Einstein tensor for g a p, T a $ is the energy momentum tensor and 
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T a = & 



1 1 

L -Gr 2 a A F(kGa 2 )g af3 , (4) 



e aP = K(g^U M U [u>0] - V 7 -g^U [aM U [rM g aP ) - \U a U p , (5) 

where a pair of indices surrounded by brackets stands for antisymmetrization, i.e. A\ a Bm — A a Bp — ApB a . 
Similarly one derives a scalar equation that can be brought into the form 

[ M (kPhf'huhv) h^^] ;/3 = kG[g a P + (1 + e-^)U a U^]f a0 , (6) 
where h a/3 = g a @ — U a U^ and n(y) is defined by 

~ M-^(m) — ^ ^ 2 F{p) — Vi (7) 

where F = dF/dfi. The scalar field a is given by 

kGa 2 = ^{U 2 h af3 4>^(t)^). (8) 
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Note that the form of the function F(/j,) [or equivalently y(/z)] is not predicted by the theory and is essentially a 
free function. In next Section we give the form used in Bekenstein's paper. As it will turn out, however, the results 
derived here are essentially independent of the exact choice of F(/j,) and quite general. On the other hand, our final 
conclusions do depend on the choice of F, since it influences the way in which observations put constraints on the 
parameters of the theory (see for example Sec. II A). 

Finally, the vector equation is derived through variation of S with respect to U a 

KU la;l3] l0 + XU a + faG^U"^^^ = 8ttG(1 - e'^g^U^f^, (9) 

where U a = g af3 Uf3 and A is a Lagrange multiplier. These 4 equations determine A and three of the components of 
U a with the fourth being determined by the normalization of the vector field 

g al3 U a U p = -1. (10) 



A. The function F 



The function F(/i) (or equivalently of y(fJ>)) is a free function since there is no theory for it. One has large freedom 
in choosing the form of F, each to be checked on implications for cosmological models, galactic rotation curves and 
constraints from measurements in the outer solar system. Bekenstein in Ref. |12| made the following choice 



which, using Eq. J7J) leads to 



3 /i 2 (/i-2) 2 



= 3 M (4 + 2/x~4^ 2 + M 3 )+41n(l- M ) 
8 fi 2 



It can be shown that the range < fi < 1 (i.e., y > 0) is relevant for quasistationary systems and 2 < fi < oo (i.e., 
y < 0) for cosmology. For this specific choice of F(/j,), one can put a lower limit on the value of the k parameter of the 
theory so that it is not in conflict with the measured motions of planets of the outer solar system (see Ref. Sec. 
IV). On the other hand, small values of n are relevant for cosmological models. Together, these constraints indicate 
a value of k around ~ 0.03. It should be stressed that it depends on the specific choice of the form of F(/j,). 

The Newtonian limit of a spherically symmetric system has been explored in Sec. IV C of Ref. where it is 
shown that for gravitational accelerations |a|/ao S> 8it 2 /k 2 the quantity y — > oo and consequently /i — > 1. As an 
arithmetic example, at earth's and Mercury's orbit \i differs from unity by about 2 • 10 -6 and 5 • 10 -8 respectively for 
the specific choice of function F and k — 0.03. Since in this study we focus on the strong gravity limit, we can 
safely take \i = 1 and, therefore, [see Eq. iJSJ] 

Strictly speaking, /i has been shown to be of order unity in the Newtonian limit but not necessary in the relativistic 
limit. However, using the analytic solutions derived in Sec. IV, we have checked that taking fj, = 1 is an excellent 
approximation also in this limit. 



III. SPHERICAL SYMMETRIC, STATIC SPACETIMES IN TEVES 



From this point on we focus on the strong gravity limit of TeVeS and explore the spacetime in the vicinity of a 
spherically symmetric mass. The isotropic form of a spherical symmetric, static metric is 

g a pdx a dx p = -e v dt 2 + e c (dr 2 + r 2 d8 2 + r 2 sin 2 8dip 2 ), (14) 

where both is, £ are only functions of r. For the static system of our case the vector field has two non-vanishing 
components 

U a - (t/*,C/ r ,0,0), (15) 
where U l and U r are functions of the radial coordinate. 
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Taking <p = <p(r), the scalar-field equation © in vacuum may be written as 

p -(^+3C)/2 



[r a e(" +3 «/V(e- c - (CT) 2 )]' = 0, (16) 

where a prime stands for ordinary derivative with respect to r. 
Equation can be integrated once to give 

^ = C r 2 [e-C-(Tp)T M 

where C is an integration constant. To determine this constant, one must consider the source of the gravitational 
field. A analysis relevant for an extended source (e.g. a star) is given in Bekenstein ^2] Sec. V, where it is shown 
that one can define the "scalar" mass m s as a (non-negative) particular integral over p and P (defined as the proper 
energy density and pressure expressed in the physical metric) of the star's matter so that the integration constant is 
given by 

C = KGm s /(4n). (18) 

The r component of the vector equation (the 9 and <fr components vanish because of the symmetry of the problem 
under consideration) can be brought into the form 

/ K{Gm.)* e-C-"a x 
U \ X+ 2n r^-{Urff)- {) - (19) 

This equation shows that there are two cases: either U r vanishes or one has a constraint on the Lagrange multiplier 
A. In the former case the vector field is aligned with the time direction, while in the latter there is a non-vanishing 
radial component of the vector field. Since the mathematical analysis of the two cases is rather different, we examine 
them separately. 



A. Case I: the vector field is aligned to the time direction 

When U r vanishes, the vector field is determined by the normalization expression (|10|) which yields 

U a = (e"^ 2 , 0,0,0). (20) 

The physical metric is then given by Eq. J5J which reduces to gtt = gtte 2 ^, gu = gue~ 2 ^. To determine g a ^ ', one needs 
to solve for v and £ and, through them, for the metric g a p and the scalar <p. To this end, the differential equations 
resulting from the tt, rr and 99 components of the metric equation must be solved. This is equivalent to the 
procedure one follows to arrive at the GR solutions. 

Since we are looking for vacuum spacetimes, the terms that include the matter energy density in Eq. (|3J) are zero. 
So, we are left with the r Q( g and Q a p terms. The r Q( a [see Eq. (@J] contains <j). a terms while the last term depends on 
the function F. In the strong acceleration limit, it is possible to show (see Sec. V in Ref. [l^) that the last term is 
completely negligible in comparison to the other terms. This is exactly the limit we are interested in here and so we 
will neglect the F term. Using Eqs. (@J, l(T7Jl and J20J), we then find 

K,Gm 2 s e~ 2 ^ 

T tt = — , ( 21 ) 



32tt 2 



(22) 



KGm 2 s e-«+^ 

T ee = — 77^, 2 — ■ ( 23 J 

32n z r z 

To proceed with the calculation of the Q a p terms [defined by Eq. (JSJ], one first needs to compute the Lagrange 
multiplier Ai (where the index I is used to show that it corresponds to the case I), from the t component of Eq. @ 
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(the other components vanish because of the symmetry of the problem and because we study the case where U r = 0). 
Using Eq. JgJJJ and that U cj)^ = and f a(i = 0, we have 

Ax = -^( T + ^ + 7 ). (24) 
Substituting Eqs. JHJ) and JT5J) in (J5J we get 

Qrr = ~ (26) 

©00 = f (rv')\ (27) 

We can now use the tt and rr components of the metric equation [Eq. ©] to derive a system of ordinary differential 
equations for £ and v. Using Eqs. (|21|l . (|22p. (|25|) and H2bj) in J3J), and after some rearrangement, one finds that 

(" + ^L + A = ^ tl iff - — — H 1 — H ) (28) 

S 4 r 4tt r 4 \8 2 4 r/ v; 

and 

(CO 2 CV C' + ^_ K (Gm s ) 2 e^) ( y /)2 
4 2 r 4tt r 4 8 1 ' 

These two equations are in principle enough to solve the metric. However, it turns out that it is useful to make use 
also of the 99 component of the metric equation 

-^— + — + ^h- = 4^ — +K —- (30) 

The study of the properties of these equations for the appropriate boundary conditions constitutes most of the rest 
of this work. 

B. Case II: the vector field has a non- vanishing r component 

In the non-aligned case, U r ^ 0, the Lagrange multiplier An is given by [see Eq. ifH^l ] 

K,(Gm s ) 2 e-^O 
11 ~ 2tt r^[e-C-(U-) 2 } 2 ' 1 ' 

The components of the vector field are connected to the functions v and £ through the normalization Eq. (|10() 

e "(f/*) 2 - e c (U r f = 1 (32) 

and the t component of the vector equation (0 (given in a compact form) yields 

m n(Gm s ) 2 e-^O 
KU *~ 2tt r*[e-<-(Ur)*]* U " °" (33) 

where the U^'^.^ term involves derivatives of the four unknown functions v, C, U r , U l . 

The last two expressions can be combined with the tt and rr components of the metric equation Q to arrive to 
a closed system of four differential equations with four unknown functions. To this end, one has to calculate the 
relevant T a p and Q a terms. Equation ||3J yields 

nGm 2 e~ 3? 

Ttt = 32tt 2 r 4 [e-C - (CA) 2 ] ' (34) 
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(l-3(£T) 2 ). (35) 



_ nGm 2 s e -(" +3 » 
Trr ~ 32tt 2 r 4 [e-C - (U r ) 2 ] 2 

For 8 tt and 8 rr we have [see Eq. |JSJ] 

e - g e *W i m>? i K{Gms? V r ^ (v+20 07) 

The task in the next subsection is to study the asymptotic behavior of the physical metric far from the source in the 
two cases I and II and derive the Post-Newtonian corrections predicted by TeVeS. 



C. Asymptotic behavior of the metric far from the source 

Far from the source (but not too far, so that the MOND corrections can be safely neglected), the metric can be 
taken to be asymptotically flat. Expanding the e^, e v to powers of r/r g (where r g is a length scale to be determined) 
we have 

= l-r g /r + a 2 (r g /r) 2 + ■■■ (38) 

and 

e< = 1 + hr g /r + b 2 {r g /r) 2 + ■ • ■ , (39) 

where the proportionality constant of the second term in the expansion (|38|l has been absorbed by r g . We now proceed 
to calculate the coefficients a, and bi of the metric and equivalent coefficients of the physical metric g a p for the two 
cases I, II (defined in the previous Section). 



1. Case I: U r = 



If the vector field is aligned with the time direction one can substitute the expansions l|38|) . 
equations (|28|) and (|29[1 . match coefficients of like powers if 1/r and solve for the coefficients aj, 
of (1/r) 3 , the metric has the form 



r 



1 r 

2r 2 



and 



1+^ + TTT 



2k / Gm s \ 2 
7r V r„ 



1 

9G 



K 



18+ — 



Gi 



K 



13911 into the metric 
bi. Doing so to order 



(40) 



1 

96 



10k fGr 



7T V 



5K 



(41) 



In this expansion, one can see that the first corrections introduced by TeVeS with respect to the Schwarzschild metric 
appear in the {r g /r) 2 term in and in the {r g /r) 3 term in e". 

Actually, these asymptotic expansions differ from expressions given in [l2j (compare Eqs. (|40|l and (|41|) of this work 
with Eqs. (89)-(91) in Section V of [l2^h The reason for this difference is a sign error in the /3i term of the Lagrange 
multiplier in Bekenstein's Eq. (82) (see also the erratum of Ref. 01 )■ Because of this discrepancy, we need to rederive 
the Post-Newtonian corrections predicted by TeVeS. The physical metric g a p is given by the expressions g tt — g tt e 



9u — gue so we still need the asymptotic behavior of i 



nGm s 
Aitr 



KGm s ' 

1927T 



1 



2<f> 

Integrating l|17l) and using Eqs. (|40|l and l|41|). we have 
Gm 



K 
~2 



(42) 



where <p c is the cosmological value of at a specific epoch, which can be absorbed by rescaling of the t and r 
coordinates: t' = te^ c and r' = re~* c . Doing so and dropping the primes for simplicity in the notation, the physical 
metric is 



gtt = -1- 



nGm s 
2tt 



1 



2 1 



nGrrit 



)K 



K,Gm s 



2k iGm s \ 2 
7r V r„ 



-K 



, r 2 

9 



(43) 



7 



and 

(kGvh s \ 1 1 r (K,Gm s \ 2 „ 2k /Gm s \ 2 „ir 

.9rr = l+Hr^ +r 9 - + T^ 2 + 2 -2-— A + #-§+••■ 44 

V Z7r / r 16 L V 7rr 9 / tt \ r g / 1 r z 

Identifying the \ jr term of the tt component with 2Gjqmlr (where Gat is Newton's constant), the physical metric 
can be brought into the form 

gu = -1 + 2G N m/r - 2G 2 N m 2 /r 2 + ■■■ (45) 

and 

g rr = 1 + 2G N m/r H (46) 

which is identical to GR up to order of post-Newtonian corrections. This means that one has to go to higher order 
terms in TeVeS to obtain the corrections to GR. No constrains can be set to the parameters of TeVeS K and K from 
measurements of the standard Post-Newtonian coefficients, if the radial component of the vector field vanishes. 

At this point, one more comment is in order. By inspection of Eqs. I|4()jl. (|41(l one notices that the quantity 
^(Gm s /r g ) 2 — K/2 (times some factor) appears in all the corrections introduced by TeVeS with respect to equivalent 
the general relativistic solution. This quantity will also appear in the analytic solutions derived in Sec. V. 

2. Case II: U r / 

In the non-aligned case one needs to consider the asymptotic expansion of the vector field components. For r g /r -C 1 
the vector field relaxes to its cosmological value, i.e., J7* — > 1 and U r — > (since there is no preferred spatial direction). 
So, expanding to powers of r g /r, we have 

U l = l + Cl r g /r + c 2 (r g /r) 2 + ■■■ (47) 

and 

U r = d irg /r + d 2 (r g /r) 2 + ■■■. (48) 

From this point on the method we follow to calculate the Post-Newtonian corrections is similar to this of the previous 
subsection. Substituting these expansions and Eqs. I|38|l . (|39|l into Eqs. (|32J) . and the tt and rr components of the 
metric equations © and matching coefficients of like powers if 1/r, we derive the coefficients dj, bi, a and di. This 
analysis is carried out down to the order necessary to calculate the post-Newtonian coefficients and gives for K <C 1, 

K < 1 

1 K<Gm s ) 2 K 



2 Airr 2 g 



1, (50) 



, 3 3 K (Gm s ) 2 K 

b2= a~a 2 77^' ( 51 ) 

8 8 nri 16 



ci = \, (52) 



C2 



1 / lK{Gm 8 y Sn{Gm s f x 

16 1 5 + Knr* ± V Knrj + V (53) 



and 



1 / 8K(Gm s ) 2 \ , , 

di = i( 1± v-bf L+5 )- (54) 
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The ± sign in the last two coefficients comes from the fact that the normalization expression ilML'l contains squares of 
the vector components. The + sign corresponds to U r > and vice versa. 

The asymptotic behavior of the scalar field is found after expanding and integrating Eq. (|17|1 

<P(r) = <Pc-^ + 0(r- 3 ). (55) 

The physical metric is given by Eq. J2J) and a rescaling of the r, t coordinates by t' = te^ c and r' — re~^ c is needed 
so that it can asymptote to the Minkowskian form. Notice, however, that </> c is not absorbed by this rescaling unlike 
the U r = case, because of the more complicated connection of the physical metric with the fields of TeVeS in this 
case. Assuming again K <C 1, K <C 1 and furthermore that 4> c <C 1, we have for the standard /3, 7 Post-Newtonian 
coefficients, as predicted by TeVeS for the case that U r 7^ 

^ 1+ 8^ + 4 + ^( 3 + ^ ± fe + 5 ) (56) 

and 



7=1. (57) 

Here again the ± sign in the expression for (3 is determined by the sign of U r [see Eq. (I54ll ]. 

While the 7 coefficient coincides with the GR prediction, the (3 differs from unity. The best determination of (3 
comes from lunar laser ranging tests (see for example Ref. |1 5| 't which in combination with the value for 7 measured 
by the Cassini experiment [Tg yields (3 — 1 < 10~ 4 How does this result compare with the prediction of Eq. H56fl ? 
It is important to note that the term multiplied with the (positive) 4> c in Eq. I156H is positive for any value of k/K > 
and choice of the ± sign, so one has the inequality — 1 > k/(8ty) + K/-L. For the choice of the function F made in 
Ref. K is constrained to be ~ 0.03 which results in (3 — 1 > 2.5 ■ 10~ 3 (taking the K term much smaller). This is 
in conflict with observations. 

Summarizing, in this section we have shown that if the vector field is aligned with the time direction, the standard 
Post-Newtonian coefficients derived by TeVeS are identical to these of GR, while if U r 7^ 0, the PPN correction for 
the j3 coefficient is in conflict with best determinations of (3. This means that either one has to assume that U r = 0, 
or a different choice of the function F than that of Ref. [IJ has to be made so that TeVeS is in accordance with solar 
system phenomenology. 



IV. ANALYTIC SOLUTIONS WHEN U r VANISHES 



Until now, we have kept the study of spherical symmetric spacetimes in TeVeS quite general. From this point on, 
we focus on the branch of solutions for which U r — 0, i.e., the vector field is aligned to the time direction. As it turns 
out, exact analytic solutions are possible in this case. 



A. Solutions in the K — * limit 



The system of Eqs. (I28L (I29|) and l|30fl is rather complicated. Here, we first consider some special cases and then 
use the intuition we gain to derive the general solution. In the simplest case where both k = and K — 0, the metric 
equations in TeVeS coincide with these in GR and their write hand side is zero (i.e. no source terms appear). In this 
limit the integration of Eqs. (|28[l . (|29|) is straightforward leading to the familiar GR solution 



V 1 + r 5 /4r/ 



e c = (l + r g /4r)\ (59) 

where the boundary conditions (|40|l . I|41() have been used. In this case, one can show [see Eq. (|17|) ] that (j> is 
constant at its cosmological value <p c and that g a p coincides with that predicted by GR. The physical metric is given 
by Qtt = 5tte 2< ^ c , gu = gue~ 2 ^ c . The factors e ±2c ^ c can be absorbed by an appropriate rescaling of the t and r 
coordinates, resulting in a physical metric is equivalent to that of GR. 
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As a next step toward the most general solution, we take the limit K — > but allow k to be arbitrary. In this limit 
we essentially decouple the theory from the vector field and the metric equations become 

4 r 47r r 4 v ' 



(C') 2 + CV + C+^ = 



4 2 r 47r 
and 

C" + (i/) 2 C' + v' _ n{Gm s ) 2 e -^ +y ) 

2 + 4 + 2r ~~ 4^ r 3 



(62) 



It turns out that that Eqs. H6U[) . (|61|) . I|62|) are equivalent to spherical symmetric spacetimes in metric-massless scalar 
theories of gravity. The exact solution was originally written down by Buchdahl in Ref. (see also Ref. ^H). Here, 
we will briefly repeat the derivation. 

From the addition of Eqs. (|6*TJ) and Ip^l we find 



which can be integrated once to give 



2(z/' + C") + (*/ + C') 2 +6^±^ = 0, (63) 



4r 

^ + C'= , 2 c 2V (64) 



Where we have introduced the integration constant r\. This constant can be evaluated by expanding Eq. (|64|) to 
powers of 1/r and comparing with the expansions (|4(J[I and i|41|) . After some algebra we find 



(65) 



Equation l|64|) can be integrated again to yield 

v + ( = 2ln[- ^), (66) 



where the second integration constant has been set to unity so that the asymptotic form of the solution is a flat 
spacetime (i.e. e y+ ^ — > 1 for r — > oc). 
One verifies that after setting 

/ 4r 2 7V 

C / 9 o\ 9 9 (67) 

r(r z — rj) r z — rj 

and using Eq. I|66l) to derive i/', the metric equations are all satisfied. After integrating for v and £, one has the exact 
solution for the metric components 



fir- — T \ r g /2r c 

\r + TV / 



and 



e C = ^^M — , (69) 



C (r 2 - r^) 2 /r - - r s/ 2r ' 
>r + r, 

where r c is given by Eq. (|65|l . It is straightforward to check that in the limit where Km s — > 0, one derives the well 
known general relativistic expressions. 

Having solved for the metric components, one can integrate Eq. (|17() to derive the r dependence of the scalar field 
and then the physical metric through Eq. (J2J. However, the results derived in this Section are of limited generality 
since they correspond to the K = case, where the effect of the vector field to the metric equations is ignored. The 
generalization of the solutions to the case where K ^ is the task of the next Section. 
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B. Spherically symmetric, vacuum solution for the metric for arbitrary k, K 



We turn to the general case where both n and K are non-zero. While at first sight the metric equations look quite 
complicated in this case, it turns out that one can repeat the procedure of the previous Section to derive more general 
spherical symmetric, vacuum solutions for the metric which are identical to l|t)8|) . I|t)9|) provided that one makes the 
substitution 



in the definition of r c [Eq. (|u5f) ]. i.e. 



In the k, K parameter plane the line defined by 

K 



2*(Gm.y 
7r V r g J 



has a special significance. With Eq. Q72JI T c — r g and the metric is identical to the general relativistic one |l8j . This 
is a result of the fact that the energy density contributed by the scalar is exactly canceled out by the negative energy 
density of the Qtt in the right hand side of the tt component of Eq. (J3J), i.e. &itGtu + Qtt = 0. Actually, when 

2k ( Gm s \ 2 . „ 

K >-{ — ) (73) 



tt v r g 

the total energy density of vacuum contributed by the fields is negative in the whole space time. This can have 
important consequences for the theory since it may lead to instability of the vacuum from the quantum point of view. 
The behavior of the scalar field can be followed by integrating Eq. 1|17[) and the use of Eq. (jfifijl 

<P(r =0 C +— -^ln — -i , (74) 
oirr c \r + r c / 

where <p c stands for the cosmological value of the scalar field at a specific epoch. Just as in metric-scalar theories, 
one can see that, unless nm s = the scalar field diverges logarithmically at r = r c and that there is always a radius 
r\ > r c where 4>(ri) — and becomes negative further in. We will return to this point in the next Section where we 
discuss how black holes look like in TeVeS, and how much they differ from the ones predicted by GR. 

The components of physical metric is related to v, Q and 4> through the expressions @, i|68|) . and I|tj9|) that yield 



9tt : 
and 

(r 2 — r 2 ) 2 fr — r, 



r-^-r (75) 

V r + r r / 



(76) 



K-Gm s 



where a = ^ , 47rrc . 

The expressions l|75|l and l|76|l describe spherically symmetric, vacuum spacetimes, i.e., they describe the spacetime 
down to the surface of a star. The two dimensionless parameters k and K of the theory provide the parameter space 
that is to be explored. In addition to these, we have kept the scalar mass m s and the gravitational radius r g as free 
parameters so that the derived results to be quite general and applicable to both the case of a black hole and the 
exterior of a star. In appendix D of Ref. [T^ a detailed description of the procedure to calculate m s and r g in terms 
of its gravitational mass m g of the star is given. Unfortunately, this method is not applicable to the case of a black 
hole and a different approach is needed to determine m s and r„. 
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V. HOW DO BLACK HOLES LOOK LIKE IN TEVES? 



The characteristic radius of the physical metric described by Eqs. (|75|l and (|76|l is r c . At r c the tt component of the 
metric vanishes and the question that rises is whether and under which conditions r c is the location of the horizon of a 
black hole. A first step toward answering this question is to calculate the surface area at this radius, which turns out 
to be proportional to g rr {i"c)- A black hole must have a finite surface area at r c which [in view of Eq. I|76(l ] constrains 
a to be < 2. 

A second constraint on a (Bekenstein private communication) comes from the demand that there is no essential 
singularity at r c . For our solution, the Ricci scalar R is 

2(a 2 -4)r c V(r-r c r- 4 
R (r + r c )«+4 ■ ( 77 ) 

From this expression one can see that the Ricci scalar is finite when a = 2 or a > 4. Considered together, the two 
constraints (i.e. of finite surface area and Ricci scalar at r c ) imply that only the value a = 2 describes a black hole. 
Using the definition of a we have that for a = 2 



nGm s 
4 ' 8vr 

and the physical metric has the form 



(78) 



~9tt = -{— ± ) (79) 
\r + r c J 

and 

ffrr=(l + ^) 4 . (80) 

This is exactly the GR solution after setting r c = GnUi/I. So, the physical metric in TeVeS is identical to this of GR 
for a non-rotating black hole. 

Furthermore, one can use the definition of r c [see Eq. I|71f) ] in Eq. (|78|l to solve for m s and finds 



1+ Wl+ (2-&)2£ 



Gm. V V 27r 



Z 2tt 



(81) 



We have already shown in the previous Section that when m s 0, there is always a region close to r c where the scalar 
field turns negative. Bekenstein in ref. 011 the other hand, has shown that TeVeS becomes acausal (i.e. it suffers 
from superluminal propagation of metric, vector and scalar field disturbances) when (f> < 0. As a result, the theory 
appears to behave in an unphysical way in the vicinity of our black hole solution. On the other hand, our solutions 
have been derived under the assumption that U r — 0. Perhaps the causality problem can be overcome by allowing 
for U r ^ 0. 



VI. CONCLUSIONS 



Bckenstein's recent relativistic gravitational theory (TeVeS) that leads to MOND in the relevant limit has been 
proposed as a modification to GR. TeVeS has several attractive features, for example it predicts the right amount 
of gravitational lensing when only the observed mass is used and provides a covariant formulation to construct 
cosmological models. 

The free parameters in TeVeS can be constrained by the large extra-galactic phenomenology. In this work, instead, 
we have looked at TeVeS in the strong gravity limit. Two branches of solutions are identified: the first is characterized 
by the vector field being aligned with the time direction while in the (not previously explored) second branch the 
vector filed has a non- vanishing radial component. We have shown that the (3 and 7 PPN coefficients in TeVeS are 
identical to these of GR in the first branch of solutions while the (3 PPN coefficient differs in the two theories in 
the second. Despite the fact that the results derived here are essentially independent of the exact choice of the free 
function F of the theory, our final conclusions do depend on it, since the choice of F influences the way in which 
observations put constraints on the parameters of the theory. For the second branch of solutions and for the choice 
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of the free function F made in Ref. [12|, TeVeS predicts (3 that is in conflict with recent observational determinations 
of it. 

For the first branch of solutions, we derive analytic expression for the physical metric. These solutions are an 
extension of these that describe spherical symmetric spacetimes in tensor-massless scalar theories and depend on the 
values of the two dimensionless parameters k, K of TeVeS and the ratio Gm s /r g . One of the findings of this work 
is that the energy density contributed by the vector field is negative, and when K > ^(Gm s /rg) 2 the total energy 
density of vacuum also becomes negative, possibly turning it unstable from the quantum point of view. 

In the case of a black hole, our solutions for the metric are identical to the Schwarzschild solution in GR. On the 
other hand, these solutions are shown to be acausal in the vicinity of the black hole. Possibly, the issues of the negative 
energy density contributed by the vector field, and of causality close to a black hole do not appear in the case where 
U r 7^ 0. In this case, however, a different choice of the free function F will be needed so that TeVeS is not in conflict 
with solar system phenomenology. 
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